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Lithium  iron  phosphate  is  a  promising  cathode  material  for  the  use  in  hybrid  electrical  vehicles  (HEV) 
meeting  the  demands  of  good  stability  during  cycling  and  safe  operation  due  to  reduced  risk  of  thermal 
runaway.  However,  slow  solid  state  diffusion  and  poor  electrical  conductivity  reduce  power  capability.  For 
further  improvement,  the  identification  of  the  rate  determining  processes  is  necessary.  Electrochemical 
impedance  spectroscopy  (EIS)  has  proven  to  be  a  powerful  tool  for  the  characterization  of  electrochemical 
systems.  In  this  contribution  a  deconvolution  of  the  impedance  with  the  distribution  of  relaxation  times 
(DRTs)  is  used  to  obtain  a  better  resolution  in  frequency  domain.  Therewith,  the  relevant  loss  processes 
are  identified  and  an  impedance  model  is  developed.  Using  DRT  and  CNLS-fit  allows  the  determination  of 
time  constants  and  polarization  resistances  of  all  relevant  loss  processes.  Furthermore,  their  temperature 
behavior  is  studied  and  a  physical  interpretation  is  provided. 

©  2010  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  automobile  industry  has  established  rigorous  standards  in 
terms  of  security  and  reliability.  The  success  of  electro  mobility 
is  strongly  dependent  on  advances  in  lithium  ion  battery  technol¬ 
ogy.  In  recent  years,  a  spectrum  of  candidate  materials  for  battery 
electrodes  and  electrolytes  has  been  investigated.  Since  cathodes 
fall  short  of  the  theoretical  capacity  compared  with  typical  anodes, 
many  works  concentrated  on  the  knowledge  and  improvement  of 
cathode  materials  [1-5]. 

After  being  proposed  in  1997  as  a  potential  candidate  for  the 
use  in  lithium  ion  batteries,  lithium  iron  phosphate  (LiFeP04)  has 
drawn  much  attention  due  to  its  good  cycling  behavior  [1,4,6]. 
It  offers  a  relatively  large  theoretical  capacity  of  170mAh/g  and 
allows  a  cell  voltage  of  3.4  V  vs.  Li,  maximizing  energy  density 
while  minimizing  side  reactions  such  as  electrolyte  decomposition 
[2,6].  Furthermore,  analysis  by  differential  scanning  calorimetry 
(DSC)  revealed  excellent  thermal  stability  in  comparison  with  other 
established  cathode  materials  like  LiCo02  and  LiMn204[7].  How¬ 
ever,  the  positive  aspects  are  counteracted  by  low  electronic  and 
ionic  conductivities  leading  to  a  poor  high  power  capability  [2]. 
Rate  performance  as  well  as  the  charge  and  discharge  capacity  of 
LiFeP04  is  directly  affected  by  the  cell  temperature.  To  improve 
the  electrical  and  electrochemical  properties  of  a  LiFeP04  cathode 
structure,  the  rate  determining  steps  need  to  be  clearly  identified 
and  thoroughly  understood  [4,8]. 
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Along  with  cyclic  voltammetry  (CV),  electrochemical  impedance 
spectroscopy  (EIS)  has  proven  to  be  a  powerful  tool  for  character¬ 
ization  of  electrochemical  systems  in  general  [1,8,9].  Impedance 
spectra,  which  contain  information  about  the  physico-chemical 
processes  within  the  cell,  are  usually  evaluated  by  an  equivalent 
circuit  model.  This  approach  requires  knowledge  about  number, 
and  physical  origin  of  all  processes,  that  contribute  to  the  cumu¬ 
lative  impedance.  Fig.  1  gives  an  qualitative  overview  of  processes 
adding  to  ohmic,  polarization  as  well  as  diffusion  resistance  in  a 
cell  with  LiFeP04  cathode  structure  and  the  order  of  magnitude 
of  their  expected  frequency  range  as  given  in  [10,11].  Processes, 
which  can  be  assigned  to  the  lithium  metal  counter  electrode, 
are  not  considered  in  this  figure.  As  these  processes  lower  the 
voltage  of  the  cell  by  causing  an  overvoltage  during  operation, 
these  processes  will  be  further  considered  as  loss  processes.  The 
identification  even  of  the  most  prominent  loss  processes  and  the 
determination  of  their  specific  time  constants  remain  ambigu¬ 
ous,  as  they  usually  overlap.  This  challenge  can  be  addressed  by 
a  deconvolution  of  the  impedance  with  the  distribution  of  relax¬ 
ation  times  (DRTs).  This  advanced  method  has  been  developed  for 
and  successfully  applied  in  high  temperature  fuel  cell  research, 
where  a  physically  motivated  impedance  model  without  the  need 
of  any  a  priori  knowledge  of  the  investigated  electrochemical  sys¬ 
tem  was  generated.  The  DRT  method  offers  a  higher  resolution  in 
the  frequency  domain,  allowing  a  clearer  identification  of  the  loss 
processes  [12-14]. 

In  the  present  work,  three  cell  configurations,  (a)  full  cells 
(LiFeP04/Li),  (b)  symmetrical  cathode  cells  (LiFeP04/LiFeP04)  and 
(c)  symmetrical  anode  cells  (Li/Li)  are  evaluated.  The  symmetri¬ 
cal  setup  allows  to  separate  anodic  and  cathodic  loss  processes.  All 
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Fig.  1.  Most  prominent  loss  processes  in  a  lithium  ion  cell  with  a  LiFeP04  cathode  structure,  sorted  by  their  potential  characteristic  frequencies. 


test  cells  are  analyzed  by  impedance  spectroscopy  in  a  tempera¬ 
ture  range  of  0-40  °C.  The  low  frequency  branch  of  the  ElS-spectra 
is  modeled  by  a  series  of  a  capacity  C  and  a  generalized  finite 
length  Warburg  (GFLW)  and  then  subtracted.  For  the  remaining 
impedance  data  the  distribution  of  relaxation  times  (DRTs)  is  calcu¬ 
lated.  This  two-step  pre-identification  of  all  relevant  loss  processes 
and  their  time  constants,  leads  to  an  equivalent  circuit  model, 
which  is  finally  used  to  derive  the  area  specific  resistance  (ASR) 
and  activation  energy  (Efl)  for  all  processes  found.  This  methodol¬ 
ogy  constitutes  a  new  approach  to  investigate  lithium  ion  battery 
electrodes,  delivering  an  accurate  separation  of  two  loss  processes 
in  the  LiFeP04  cathode. 


2.  Experimental 

For  this  study,  cathodes  were  prepared  from  basic  carbon  coated 
LiFeP04  powder  (Sued-Chemie  AG),  to  which  carbon  black  and 
PVDF-binder  were  added  to  receive  a  final  weight  ratio  of  70:24:6 
(LiFeP04:carbon  blaclcbinder).  The  mass  of  active  material  in  each 
cathode  was  about  1.5  mg.  Thereafter,  a  l-methyl-2-pyrrolidinone 
(NMP)  based  slurry  was  prepared,  applied  on  an  aluminum  foil  by 
doctor-coating  and  vacuum  dried  for  120  min  at  80  °C.  The  elec¬ 
trodes  were  punched  out  and  finally  tempered  for  another  180  min 
at  120  °C. 

All  measurements  were  conducted  on  Swagelok-type  test  cells 
in  three  different  configurations:  full  cells  with  LiFeP04  vs.  metal¬ 
lic  lithium  (a),  as  well  as  symmetrical  setups  with  LiFeP04  vs. 
LiFeP04  (b)  and  lithium  vs.  lithium  (c).  All  investigated  cells  fea¬ 
tured  electrodes  with  1.2  cm  in  diameter  and  an  active  area  of 
1.13  cm2  respectively.  The  metallic  lithium  foil  (Sigma-Aldrich) 
had  a  thickness  of  0.38  mm  which  resulted  in  a  mass  of  active  mate¬ 
rial  around  23.3  mg.  This  guaranteed  a  sufficiently  high  lithium 
reservoir  for  the  examination  of  cathode  material.  The  separators 
of  all  cells  consisted  of  two  layers,  a  Freudenberg  FS201 9  serving  as 
an  electrolyte  reservoir  and  a  Celgard  C500  membrane,  to  prevent 
short-circuiting.  The  electrolyte  used  was  a  1  mol  lithium  perchlo¬ 
rate  (LiC104)  solution,  based  on  ethylene  carbonate  (EC)  and  ethyl 
methyl  carbonate  (EMC)  with  a  gravimetric  mixing  ratio  of  1 :1. 

After  cell  assembly,  ten  charge  and  discharge  cycles  were  per¬ 
formed  to  ensure  a  steady  state  of  the  system  (compare  Fig.  3). 
Impedance  measurements  were  carried  out  using  the  Solartron 
1 400E  cell  test  system  with  Scribner  Multistat  software.  The  applied 
ac  voltage  was  10  mV  (RMS)  at  open  circuit  condition  and  the 
frequency  was  varied  within  a  range  of  1  MHz-10  mHz.  The  qual¬ 
ity  of  the  obtained  impedance  data  was  analyzed  by  calculating 
the  Kramers-Kronig  residuals.  Within  the  frequency  range  from 
1 0  mHz  up  to  1 00  kHz,  the  residuals  (real  and  imaginary  parts)  were 
below  0.5%,  proving  high  data  quality  of  all  conducted  impedance 
measurements.  Therefore  only  the  frequency  range  from  1 00  kHz  to 
10  mHz  was  used  for  further  analysis.  All  impedance  spectra  were 
measured  at  a  state  of  charge  of  100%  for  the  cell  system,  which 
means  a  fully  deintercalated  state  for  the  LiFeP04,  varying  temper¬ 
ature  between  0 0  C  and  40 0  C  in  a  climate  test  chamber.  At  least 
four  identical  cells  have  been  measured  in  parallel  to  validate  the 
reproducibility  of  the  results. 


3.  DRT  and  preprocessing 


The  distribution  of  relaxation  times  (DRTs)  is  a  recommendable 
method  to  deconvolute  impedance  data,  because  of  its  outstanding 
potential  to  resolve  polarization  processes  with  close-up  time  con¬ 
stants  [12-14].  Contrary  to  a  fit  procedure  with  equivalent  circuits, 
the  individual  processes  that  contribute  to  the  overall  impedance 
are  represented  without  any  assumption  for  their  physical  origins. 
Although  the  DRT  has  proven  to  be  a  useful  tool  for  the  investiga¬ 
tion  of  high  temperature  fuel  cells  [12-14]  and  has  been  examined 
by  other  researchers  [15-17]  it  remained  widely  unknown  -  prob¬ 
ably  because  of  its  high  portion  of  signal  processing.  The  relation 
between  the  impedance  and  its  DRT  is  given  by 

Z(co)  =  R0+Zpol(co)  =  R0+Rpol (1) 


with  the  condition  that 


g(r)dr  =  1 


(2) 


where  Z((o)  is  the  impedance  data,  R0  is  the  ohmic  part  of  the 
impedance,  Zpoi{co)  is  the  polarization  part,  Rpo/  is  the  polarization 
resistance  of  the  impedance  andg(r)  is  the  corresponding  distribu¬ 
tion  of  relaxation  times.  The  expression  in  the  integral  resembles 
the  definition  of  an  RC-element: 


Zrc{u>)  = 


R 

1  +j(JL>  X 


R 

1  +jooRC 


(3) 


If  the  distribution  of  relaxation  times  g(  r)  is  a  sum  of  Dirac-pulses 


N 

g(r)  =  y^<S(F  -  T„) 
n- 1 


(4) 


the  integral  represents  a  series  connection  of  N  RC-elements  with 
individual  relaxation  time  constants  r2.  Fig.  2a  demonstrates  the 
advantage  of  the  DRT  for  the  simple  case  of  N  =  2.  In  the  Nyquist  plot 
the  two  RC-arcs  hardly  can  be  distinguished  while  its  representa¬ 
tion  in  the  DRT  clearly  shows  two  Dirac  impulses  corresponding 
to  two  RC-elements.  The  height  of  the  peak  equals  the  polariza¬ 
tion  R  and  its  position  sets  the  relaxation  time  r2  or  characteristic 
frequency. 

RC-elements  describe  an  “ideal”  system  with  lumped  param¬ 
eters  whereas  RQ-elements  are  an  adequate  representation  for 
“real”  electrochemical  systems  with  distributed  parameters.  The 
RQ-element  is  here  defined  as: 

Zro(co)  = - - - jr  (5) 

u  1+ORQf 

Fig.  2b  shows  the  DRT  of  2  RQ-elements  with  a  broadened  distri¬ 
bution  curve  and  two  maxima  on  mean  value  of  the  distribution 
of  relaxation  times.  The  integral  over  a  single  peak  represents  the 
polarization  R  in  the  definition  of  the  RQ-element,  its  center  the 
relaxation  time 


(6) 
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Fig.  2.  Nyquist  plot  (center)  and  distribution  function  of  relaxation  times  (right)  for  two  different  systems:  (a)  lumped  parameter  model  with  two  series  RC-elements,  (b) 
distributed  parameter  model  with  two  series  RQ-elements. 


The  DRT  representation  of  the  2  RQ-elements  gives  a  clear  proof 
for  its  superiority,  as  both  processes  are  separated  and  their  char¬ 
acteristic  frequencies  are  apparent.  For  the  first  time,  the  DRT  shall 
now  proof  its  usefulness  to  the  separation  of  electrochemical  loss 
processes  in  lithium  ion  cells.  However,  in  contrast  to  high  temper¬ 
ature  fuel  cells,  an  extension  of  the  DRT-method  is  essential  for  the 
investigation  of  lithium  ion  cells. 

The  distribution  of  relaxation  times  g(  r )  in  Eq.  ( 1 )  is  expressed  as 
an  integral  from  zero  to  infinity.  Furthermore,  impedance  measure¬ 
ments  of  lithium  ion  cells  never  cover  the  entire  frequency  range. 
This  problem  dissolves  if  all  processes  lie  inside  the  measured 
frequency  range.  However,  the  low  frequency  processes  in  the 
mHz-range  (e.g.  solid  state  diffusion)  are  (at  least  partly)  excluded 
due  to  the  long  measurement  time. 

Another  difficulty  is  that  the  impedance  spectra  of  lithium  ion 
cells  always  include  a  series  capacity  or  a  differential  intercala¬ 
tion  capacity,  as  explained  in  [10].  As  the  distribution  of  relaxation 
times  g(r)  in  Eq.  1  does  not  include  pure  capacitive  behavior,  it  is 
mandatory  to  eliminate  the  series  capacity  in  a  preprocessing  step. 

The  differential  intercalation  capacity  and  the  solid  state  dif¬ 
fusion  are  divided  from  the  entire  impedance  expression  by  the 


Cycle  Number 


following  procedure,  consisting  of  three  consecutive  steps: 

1.  Modeling  the  low  frequency  branch  (preprocessing). 

2.  Subtraction  of  the  low  frequency  branch  (preprocessing). 

3.  Evaluation  of  the  resulting  spectra  by  DRT-approach  (main  step). 


As  a  prerequisite,  the  low  frequency  branch  must  be  physically 
reasonable  modeled,  fitted  and  subtracted.  In  agreement  with  lit¬ 
erature,  we  assigned  the  low  frequency  branch  to  (i)  solid  state 
diffusion  in  the  intercalation  electrode  and  (ii)  to  intercalation 
capacity  [10,18,19].  Modeling  approaches  are  described  in  [10,20], 
we  made  use  of  the  diffusion  and  intercalation  model  proposed  by 
Levi  and  Aurbach  [10]  comprising  of  a  finite  length  Warburg  ele¬ 
ment  and  a  series  capacity.  This  model  approach  was  also  justified 
by  the  theoretical  work  of  Jacobsen  and  West  [21]  for  diffusion  in 
a  spherical  particle  with  diffusion  towards  the  center.  Herein  the 
behaviour  at  high  frequencies  is  equivalent  to  a  finite  length  War¬ 
burg  impedance  in  series  with  a  capacity.  This  modeling  approach 
leads  to  the  impedance 


ZDiff  =ZW  +ZC  =  Rw  ■ 


tanh([j'r1A/o;|Pw)  1 
|jTw®fw  +  jcoC0 


(7) 


Fig.  4a  shows  step  1  of  the  preprocessing:  the  resulting  CNLS 
fit  of  this  model  to  the  measured  impedance  curve.  Accordingly, 
Fig.  4b  shows  step  2  of  the  preprocessing:  it  is  now  subtracted  from 
the  measurement  data  and  brings  out  an  impedance  that  can  be 
deconvoluted.  Implicitly,  two  equations  have  been  applied: 


1.  The  (total)  impedance  is  assumed  as  one  part  that  can  be 
described  by  the  DRT  and  another  part  resulting  from  solid  state 
diffusion/intercalation  and  cell  capacity 


ZtotO'®)  =  R*0, 


g(r) 

1  +  joox 


dr  +  Zflw  +  Zc 


(8) 


2.  Solid  state  diffusion/intercalation  capacity  is  subtracted  from  the 
measured  (total)  impedance 


ZpreijO)) 


=  Ztot  -  ZFlw  -  Zc  =  R*ol 


g{r) 

1  +j(JL)  X 


dr 


(9) 


4.  Results  and  discussion 


Fig.  3.  Specific  discharge  capacity  of  one  LiFeP04  cell  during  the  formation,  corre-  . 

sponding  to  its  mass  of  active  material.  The  capacity  is  very  stable  which  indicates  Fig.  5a  shows  on  the  left  hand  Side  the  preprocessed  impedance 

a  steady  state.  curves  of  LiFeP04/Li-cells  for  varying  temperatures  between  0 0  C 
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Fig.  4.  Different  stages  of  preprocessing  of  impedance  data:  (a)  fitting  of  diffusion 
and  intercalation  model  to  the  impedance  data  and  subsequent  subtraction;  (b) 
preprocessed  impedance  to  be  evaluated  by  the  DRT. 


and  40  °C.  The  impedance  curves  reveal  an  increasing  overall 
polarization  and  ohmic  resistance  for  decreasing  temperatures. 
Furthermore,  they  show  an  increasing  number  of  polarization  pro¬ 
cesses.  For  the  0 0  C  measurement,  three  polarization  processes  can 
be  identified.  Fig.  5a  shows  on  the  right  the  corresponding  DRT  of 
these  impedance  curves.  Even  for  40 0  C  it  is  possible  to  separate  five 
polarization  processes  by  the  DRT,  whereas  the  impedance  curve 
only  reveals  two  distinctive  processes. 

Over  the  entire  temperature  range,  there  are  three  major  loss 
processes  which  are  examined  in  the  following  study.  The  DRT 
illustrates  that  process  P\C  has  an  increasing  polarization  and 
a  decreasing  characteristic  frequency  (10-0.2  FIz)  for  decreasing 
temperatures.  This  behavior  is  also  shown  for  process  P^A,  whereas 
this  process  plays  a  major  role  in  the  overall  polarization.  Process 
P2c  at  1  kFIz  represents  the  main  loss  process  for  most  temper¬ 
atures  and  is  almost  temperature  independent.  The  two  small 
high  frequency  processes  cause  a  minor  part  of  the  overall  cell 
impedance  and  their  frequencies  overlap  for  low  temperatures.  For 
a  consistent  evaluation  of  these  two  processes,  a  detailed  exami¬ 
nation  in  further  studies  is  required.  In  this  first  approach  they  are 
neglected. 

For  an  assignment  of  these  three  major  loss  processes  to  cath¬ 
ode  and  anode  sides,  symmetrical  test  cells  were  assembled.  The 
impedance  measurements  of  the  symmetrical  LiFePOq  cell  are 
shown  in  Fig.  5b.  The  temperature  variation  exhibits  two  major  loss 
processes.  The  DRT  in  Fig.  5b  on  the  right  hand  side  deconvolutes 
very  clearly  these  two  loss  processes  and  their  temperature  depen¬ 
dency.  The  low  frequency  process  contributes  the  minor  part  to 
the  polarization  of  the  cathode  and  shows  a  strongly  varying  char¬ 
acteristic  frequency  (0.2-2  FIz).  The  other  process  has  a  relaxation 
frequency  around  2  kFIz  and  is  almost  temperature  independent. 

The  impedance  curves  for  the  symmetrical  Li/Li-cells  are  shown 
in  Fig.  5c  on  the  left  hand  side.  The  polarization  of  the  impedance 
curve  is  increasing  significantly  with  decreasing  temperature  and 
one  major  loss  process  can  be  identified.  The  corresponding  DRT  in 
figure  Fig.  5c  on  the  right  shows  that  the  characteristic  frequency  of 
this  process  is  shifting  from  800  FIz  down  to  100  FIz.  Furthermore, 
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Fig.  5.  Pre-processed  impedance  curves  (left  column)  for  the  three  investigated  cell  configurations:  (a)  LiFeP04/Li-cells,  (b)  symmetrical  LiFeP04/LiFeP04-cells  and  (c) 
symmetrical  lithium/lithium-cells  and  their  corresponding  distributions  of  relaxation  times  (right  column). 
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for  low  frequencies,  the  DRT  shows  an  indefinite  number  of  loss 
processes.  In  the  impedance  curve  this  is  represented  by  the  flat 
low  frequency  part  of  the  semi  circle.  This  low  frequency  domain 
must  be  examined  in  detail  to  understand  the  Li-anode  behavior. 
However,  in  this  study,  the  Li/Li-cell  is  used  for  the  separation  of 
anodic  and  cathodic  loss  processes  which  is  done  in  the  following 
section. 

4.1.  Identification  of  loss  processes  and  assignment  to  the 
electrodes 

For  this  purpose,  symmetrical  cells  and  LiFeP04/Li-cells  are 
compared  in  Fig.  5.  On  the  left  hand  side,  the  impedance  curves 
are  aligned  vertically.  The  symmetrical  cells  show  different  polar¬ 
ization  amplitudes  than  the  half  cell  impedance  curves  and  do  not 
allow  an  assignment  of  cathode  or  anode  processes.  The  right  col¬ 
umn  of  Fig.  5  shows  the  DRTs  of  all  measurements  and  illustrates 
one  main  advantage  of  the  DRT  visualization.  Processes  one  and 
three  are  directly  re-detected  in  the  symmetrical  cathode  cell  by 
comparing  their  characteristic  frequency  and  their  temperature 
behavior.  Similarly,  process  two  is  re-detected  in  the  symmetri¬ 
cal  anode  cell.  An  assignment  of  these  major  loss  processes  of  the 
LiFeP04 /Li-cell  to  anode  and  cathode  by  DRT  is  possible.  In  the  fol¬ 
lowing  sections,  processes  one,  two  and  three  are  therefore  called 
Pig  Pi  a  and  P2o 

The  exploitation  of  the  DRT’s  advantages  and  symmetrical  cell 
measurements  enables  a  clear  assignment  of  cathode  and  anode 
processes.  However,  as  mentioned  before,  the  same  processes  in 
different  cell  configuration  have  slightly  different  characteristic  fre¬ 
quencies  and  show  a  little  different  temperature  behavior.  This 
might  be  a  hint  on  difficulties  of  homogeneous  composite  layer 
fabrication  or  different  aging  statuses  of  the  measured  cells.  Nev¬ 
ertheless,  this  shows  the  manifold  prospects  of  the  impedance 
analysis  by  using  the  DRT  approach. 

4.2.  Equivalent  circuit  model 

The  preidentification  of  all  relevant  loss  processes  and  the  deter¬ 
mination  of  their  electrode  origin  lead  to  an  equivalent  circuit 
model  which  can  be  seen  in  Fig.  6.  It  consists  of  the  predetermined 
finite  length  Warburg  element  with  a  serial  capacity,  represent¬ 
ing  the  solid  state  diffusion  and  the  differential  capacity  of  the  cell 
[10].  Furthermore,  there  are  three  RQ-elements,  representing  the 
major  loss  processes  P^,  Pic  and  P2c ,  which  were  identified  by  the 
DRT  approach.  As  the  other  small  processes  P4  and  P5  only  cause 
a  minor  part  of  the  overall  cell  impedance,  they  are  neglected  in 
this  first  approach.  Lastly,  the  model  contains  a  serial  resistance 
for  the  limited  electronic  and  ionic  conductivities  of  all  cell  com¬ 
ponents  and  an  inductor  for  modeling  the  cable  inductivity.  For 
frequencies  up  to  100  kHz,  this  is  not  relevant  for  the  measured 
test  cells.  In  the  following,  this  model  is  used  for  analyzing  the 


Table  1 

Characteristic  frequency  and  area  specific  resistance  ( ASR)  of  all  identified  processes 
at  0  and  40  °C. 


Process 

fr  (Hz) 

ASR (Q  cm2) 

T  (°C) 

Ro 

- 

52 

0 

17 

40 

Pu 

30 

239 

0 

200 

19 

40 

Pdiff,C 

0.001 

1037 

0 

0.02 

110 

40 

Pic 

0.3 

50 
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Fig.  7.  Arrhenius  behavior  of  all  observed  processes,  indicating  a  reasonable  fit  result 
over  the  entire  temperature  range  from  0  to  40  °C. 

temperature  behavior  of  the  relevant  loss  processes  by  CNLS-fit  in 
Matlab®.  The  applied  starting  values  for  the  fitting  procedure  were 
directly  obtained  from  the  DRT.  The  model  was  fitted  to  impedance 
curves  for  different  temperatures  and  the  fit  results  were  evaluated. 
Table  1  contains  a  listing  of  the  temperature  specific  characteristics 
of  each  process.  Column  two  shows  the  characteristic  frequencies 
obtained  from  DRT  while  column  three  contains  the  area  specific 
resistances  (ASRs)  obtained  from  model  fit.  For  further  analysis  of 
the  obtained  model  fit  and  a  detailed  examination  of  each  process, 
the  activation  energies  are  calculated  in  the  next  section. 

4.3.  Activation  energies 

There  are  several  reasons  for  calculating  the  activation  ener¬ 
gies  of  electrochemical  processes.  The  first  one  is  checking  if  the 
obtained  fit  parameter  exhibits  Arrhenius  type  behavior  which 
indicates  a  reasonable  fit  result.  Fig.  7  provides  the  Arrhenius  plots 
of  all  polarization  processes.  The  resistance  values  follow  nicely  to 
Arrhenius  behavior  which  shows  that  the  model  fit  for  varying  tem- 


Fig.  6.  Equivalent  circuit  model  derived  from  combined  analysis  of  impedance  spectra  and  symmetrical  cell  configurations  by  using  the  DRT  method. 
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Table  2 

Comparison  of  the  determined  activation  energies  and  area  specific  resistances  (ASRs)  at  25  °  C  of  all  loss  processes  and  reported  literature  values.  Physical  origins  of  all  loss 
processes  are  given  in  the  right  column. 


Process 

ASR(25°C)  (£2  cm2) 

Pact  (eV) 

Pact, lit  (eV) 

Literature 

Physical  origin 

Ro 

28 

0.20 

0.18-0.3 

0.017 

Dasgupta  et  al.  [22] 
Chagnes  et  al.  [23] 

Electronic  conductivity  (carbon  black) 
Ionic  conductivity  (LiC104) 

Pm 

37.4 

0.47 

0.49-0.83 

Zaban  et  al.  [24] 

Charge  transfer  Li/SEI 

Pdiff,C 

207.6 

0.66 

0.4-0.95 

Takahashi  et  al.  [7] 

Maier  and  Amin  [25] 

Solid  state  diffusion  in  LiFeP04 

Pic 

16.5 

0.31 

- 

- 

Charge  transfer  between  electrolyte 
and  LiFeP04 

P2C 

176 

0.07 

- 

Gaberscek  et  al.  [8] 

Cathode/aluminum  interface 

solid  state  charge  transfer?  (P1C)  transfer  cathode/aluminum  electronic  conductivity 

diffusion  LiFeP04  /  electrolyte  through  SEI  interface  (P2C)  and  lithium  ion 

migration  in  electrolyte 
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Fig.  8.  Observed  frequency  ranges  of  the  identified  loss  processes  in  our  LiFeP04/Li-cells. 


peratures  works  consistently.  Table  2  lists  the  calculated  values  for 
activation  energies. 

The  next  step  for  exploiting  the  information  content  of  activa¬ 
tion  energies  is  the  comparison  with  literature  values.  In  this  way, 
the  reasonability  of  the  determined  values  and  the  physical  origin  of 
the  respective  process  can  be  examined.  Table  2  compares  the  cal¬ 
culated  activation  energies  for  each  process  to  reported  literature 
values. 

The  ohmic  resistance  has  several  physical  reasons.  Generally,  it 
is  caused  by  the  limited  electronic  and  ionic  conductivities  of  all 
cell  components.  Comparing  the  activation  energy  of  the  ohmic 
resistance  to  the  activation  energy  of  carbon  black  conductivity 
and  LiC104-conductivity  (see  Table  2)  gives  further  information.  The 
value  of  0.2  eV  lies  in  the  range  of  activation  energies  of  amorphous 
carbon  conductivity  [22].  This  suggests  that  the  major  part  of  the 
ohmic  losses  in  our  cathode  is  caused  by  this  low  conductivity  of 
the  cathode  composite.  However  one  has  to  keep  in  mind,  that  this 
strongly  depends  on  the  used  materials  and  the  preparation  process 
and  might  be  completely  different  for  other  LiFeP04  cells. 

The  anode  loss  process  activation  energy  being  0.47  eV  is  very 
close  to  the  values  reported  by  Zaban  et  al.  [24],  which  are  charted 
in  Table  2.  According  to  [24,26],  P^A  is  mainly  caused  by  lithium 
migration  through  several  SEI  layers. 

For  the  LiFeP04-cathode,  there  are  three  relevant  polarization 
processes  detected.  The  first  process  PdifftC  was  pre-determined  to 
be  the  solid  state  diffusion.  In  order  to  obtain  the  activation  energy 
of  solid  state  diffusion,  diffusion  coefficients  were  determined  and 
fitted  by  an  Arrhenius  equation.  A  value  of  0.66  eV  was  obtained  for 
activation  energy  which  is  between  the  reported  literature  values 
(see  Table  2). 

There  are  two  further  cathode  loss  processes  which  must  be 
examined.  The  activation  energy  for  Pic  in  the  low  frequency  region 
was  determined  to  be  0.31  eV  while  for  P2c  the  value  of  0.07  eV  was 
obtained.  In  literature,  a  value  of  0.1 5  eV  was  reported  by  Takahashi 
et  al.  [7]  for  the  high  frequency  loss  process  which  corresponds 
to  P2c  in  this  study.  However,  Takahashi  et  al.  did  not  observe  a 
second  loss  process  which  indicates  that  Pic  and  P2C  overlapped  in 
his  evaluation. 

The  identification  of  two  cathode  loss  processes  was  reported 
in  [8]  by  Gaberscek  et  al.  Gaberscek  determined  P2C  to  be  the  inter¬ 
face  between  electrode  composite  and  aluminum  current  collector 
whereas  for  Pic,  no  physical  interpretation  is  provided.  The  low 


activation  energy  being  0.07  eV  in  this  study  supports  this  physical 
interpretation  for  P2C. 

A  possible  interpretation  for  PiC  regarding  the  higher  activa¬ 
tion  energy  is  the  charge  transfer  between  electrolyte  and  cathode 
composite.  Further  investigations  are  necessary  to  confirm  this 
assumption.  A  line  up  of  all  determined  loss  processes  according 
to  their  characteristic  frequency  is  given  in  Fig.  8. 

5.  Conclusions 

Lithium  iron  phosphate  (LiFeP04)  has  drawn  much  attention  as 
cathode  material,  but  low  electronic  and  ionic  conductivities  lead 
to  a  poor  high  power  capability.  Moreover,  the  rate  determining 
steps  need  to  be  clearly  identified  and  thoroughly  understood. 

Electrochemical  impedance  spectroscopy  was  conducted  on 
Swagelok-type  test  cells  in  three  different  configurations:  full  cells 
with  LiFeP04  vs.  metallic  lithium,  symmetrical  setups  with  LiFeP04 
vs.  LiFeP04  and  lithium  vs.  lithium.  For  the  first  time,  the  evalua¬ 
tion  of  impedance  data  was  successfully  supported  by  the  DRTs 
(distribution  of  relaxation  times)  method.  This  approach  included 
a  preprocessing  by  modeling  the  low  frequency  branch  (<0.1  Hz) 
by  a  series  connection  of  capacity  (C)  and  generalized  finite  length 
Warburg  (GFLW). 

For  the  full  cell  with  LiFeP04,  we  propose  an  equivalent  circuit 
model  with  5  elements  in  series,  namely  (i)  an  ohmic  resistance  for 
electrolyte  conductivity  and  electronic  losses,  three  RQ-elements 
which  represent  the  fast  cathode  process  P2C  (cathode/aluminum 
interface  according  to  [8])  (ii),  the  charge  transfer  process  at  the 
anode  P\A  (iii)  and  the  slower  cathode  process  Pic  (charge  transfer, 
to  be  confirmed  by  further  experiments)  (iv).  The  low  frequency 
branch  Pdifjc  is  modeled  by  a  combination  of  a  finite  length  Warburg 
and  a  capacity  (v).  This  model  was  capable  of  analyzing  the  area 
specific  resitstance  (ASR)  and  activation  energy  (Fact)  separately 
for  electrolyte,  anode  and  cathode  by  CNLS-fit  in  Matlab®. 

Three  major  loss  processes  of  our  LiFeP04  cathode  structures 
were  attributed  to  physical  processes  according  to  the  indication 
attained  from  the  comparison  of  the  resulting  activation  energies 
with  literature  data: 

•  Solid  state  diffusion/intercalation:  ASR  (25  °C)  =  207.6  £2  cm2, 
E  act  =  0.66  eV. 
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•  Charge  transfer  at  the  cathode/aluminum  interface:  ASR 
(25  °C)  =  176  £2  cm2,  Eact  =  0.07  eV. 

•Pic  (probably  the  charge  transfer  at  the  cathode/electrolyte  inter¬ 
face):  ASR  (25  °C)  =  1 6.5  ^  cm2,  Eact  =  0.31  eV 
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